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Abstract 

On the basis of Quantum-Monte-Carlo results the evolution of the spectral 
weight A{k,uj) of the two-dimensional Hubbard model is studied from insu- 
lating to metallic behavior. As observed in recent photoemission experiments 
for cuprates, the electronic excitations display essentially doping- independent 
features: a quasiparticle-like dispersive narrow band of width of the order of 
the exchange interaction J and a broad valence- and conduction-band back- 
ground. The continuous evolution is traced back to one and the same many- 
body origin: the doping-dependent antiferromagnetic spin-spin correlation. 
PACS numbers: 71.30.-Fh, 74.72.-h, 79.60.-i 
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Recent results of angle resolved photoemission spectroscopy ("ARPES") revealed 
strong similarities in the low-energy excitations of a proto-type insulating copper oxide, 
i.e. Sr2Cu02Cl2 [0, and metallic cuprates like Bi 2212, Bi 2201, etc.: in both cases the 
quasiparticle (QP) band has rather small dispersion of typically 1 eV width, it is separated 
from a broad main valence-band "background" (of about 6 eV width) in much the same 
way, the /c-dispersion is similar and also the intensity modulation as function of energy is 
comparable. Thus, it appears that the dispersive metallic band evolves continuously from 
the insulating limit and has a similar physical origin as the undoped valence band in the 
cuprates. This has important consequences for the copper oxides: the excitation spectrum in 
the insulating case is decisively determined by many-body effects, documented by the known 
difficulties of one-electron bandstructure calculations for the insulating limit |l| , which then 
strongly emphasizes a many-body origin of the QP dispersion also in the metallic case. 

In this work Quantum-Monte-Carlo (QMC) results for the angular resolved photoemis- 
sion spectral weight A{k,uj) for the two-dimensional (2-D) Hubbard model are reported 
which demonstrate for this "generic" model the above strong similarities between undoped 
insulating and doped metallic situations: in particular, we find in both cases a very similar 
small dispersive low-energy band, for which the band width is set by the exchange interaction 
J ~ At'^/U when the Coulomb correlation U is of order of the non- interacting bandwidth 8t. 
This new feature is shown in the metallic case to be essentially unrenormalized (for electron- 
ically filled, i.e. u < fi states) from the insulating band, while above the chemical potential n 
additional states with again an energy scale set by J are filled in. A{k, u) is inferred from the 
QMC data by applying Bayesian probability theory in the frame of "quantified maximum 
entropy" 0,^. A consistent treatment of hyper-parameters and error-covariances along with 
high-quality QMC-data allows to reveal details of the low-lying excitations which have not 
been seen before in QMC simulations. On the basis of these QMC results and the strong 
similarities to ARPES data in both insulators and metals, we argue that the continuous evo- 
lution of the QP dispersion relation observed in the high-T^ cuprates arises from a common 
many-body origin, namely from the continuous change of the spin-spin correlation length ^ 
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and the related changes in the hole-spin correlations |p. Similarities to a phenomenology 
by Kampf and Schrieffer 0, where the magnetic correlation length is an input parameter, 
are pointed out. 

In the insulating and metallic case, respectively, the QP-like band is separated similarly 
from a broad "incoherent" background, extending about 6t to 8t. The latter corresponds 
to upper and lower Hubbard bands and follows in intensity and energy spread roughly 
the weak-coupling spin-density-wave (SDW) prediction, as already reported in earlier QMC 
simulations by Bulut, Scalapino and White 

Our results confirm, but also extend these simulation results: As stated in this previous 
work the available resolution was such that when k moved below the Fermi surface only the 6t 
to 8t broad lower Hubbard band could be resolved both for the insulating and the metallic 
P] cases. However, it was already suspected there that a narrow QP band exists, since such 



a band of width J has been found for one-hole doped in an AF insulating t-J model [|TO|,|TT 
Narrow QP bands were also observed for finite dopings in exact-diagonalization calculations 
of the one-band Hubbard |T^-|T3| and t-J models, as well as in QMC-results for the 



three-band Hubbard model P^. 



The single-band 2-D Hubbard model has the Hamiltonian 

H = -t J2 + h.c.) +UJ2 ^una (1) 

<i,j>,cr i 

on a square lattice, where t is the near-neighbor hopping, Ci^„ destroys an electron of spin a 
on site i and o- = cJ^Cj^o-. The chemical potential fi sets the filling < n >=< rij^ -|- rij j >. 
Here we focus on results obtained both for half-filling < n >=1 and other fillings in order to 
develop a systematic picture of the low-lying electronic excitations as a function of doping. 
The Coulomb correlation U is chosen equal to the bandwidth 8t and as f/ = 12t. 

In order to obtain from the QMC data for the single-particle finite-temperature propa- 
gator G{k,T) the corresponding spectral weight A{k,uj) for real frequencies u, the following 
Laplace transform has to be inverted: 



G{k,r) = -I A{k,u;)^^^^du (2) 



It is by now well established ||T8|-p^ that the maximum-entropy method provides a 
controlled way to infer the most reliable A{k, uj) in the light of the QMC data. To achieve 
the desired resolution, it is important to use a likelihood function which takes the error- 
covariance matrix of the QMC data and its statistical inaccuracy consistently into account 
The results presented in this letter are based on QMC data with good statistics, i.e. 



averages over 10^ updates of all the Hubbard-Stratanovich variables result in G{k, r)'s with 
statistical error less than 0.5%. Correlations of the data in imaginary time were considered 



by making use of the covariance matrix in the MaxEnt-procedure ||2^. As suggested in 



previous work by White [T^, various moments of the spectral weight were also incorporated 
in extracting A{k,uj). In order to check on this analytical-continuation procedure detailed 
comparisons with 4x4 exact diagonalizations |T3HT5| have been performed. 

The results presented here are for lattices 8 x 8 in size and for temperatures ranging 
from (3t = 3 (T = 0.33t) to pt = 10 (T = O.lt). Covering this temperature range allows 
us in effect to study (at half-filling) a situation where the spin-spin correlation length ^ is 
larger (for (3t = 10) than the lattice size. In this case the system behaves as if it were at 
T = and develops an AF gap. For (3t = 3, on the other hand, the spin-spin correlation 
length is shorter than our finite lattice and, consequently, the gap is diminished and metallic 
fluctuations exist As we will see, the latter situation is especially useful in interpreting 



the QMC data in the metallic (doped) regime and to relate them to those of the insulating 
(half-filled) case. 

We start in Fig. 1 by examining the single-particle spectral weight for U = 8t, (31 = 10, 
at half-filling, i.e. < n >= 1. Fig. 1(a) gives a 3-D plot of A{k, uj) versus u for fc-values out of 
the Brillouin zone, whereas Fig. 1(b) summarizes these results in the usual "bandstructure" 
u versus k plot. Here dark (white) areas correspond to a large (small) spectral weight. 

We observe in all spectra (also in the U = 12t, half-filled case, presented in Fig. 3) two 
general features: One is that A{k, uj) contains a rather dispersion-less "incoherent back- 
ground" , extending both for electronically occupied {u < ji) states and unoccupied {u > jj) 
states over ~ 6t (~ 6 eV) in the U = 8t case. The new structure, which was not previously 



resolved in QMC work is a dispersing structure at low energies with small width of the order 
of J, which defines the gap A and which, (at least for U = 12t), is well separated from the 
higher energy background. 

The splitting in the low-energy "band" and the higher-energy "background" is especially 
pronounced near r-(for u < fi) and M-(for uj > ^) points due to a relative weight shift 
from negative to positive energies as k moves through X or equally through (7r/2, 7r/2), the 
midpoint between V and M. The overall weight distribution in Aik^oj) follows roughly the 
SDW prediction as found in the QMC calculation by Bulut et al. 0: the total integrated 
weight in the SDW approximation < n^^^ >= ^'^^ A^^^ (k^Ljj)dLjj, is in good accord with 
the QMC momentum distribution 0. 

However, the dispersion of the structure near the gap does not follow the SDW prediction: 
its dispersion has a significant (about a factor of 2 for [/ = St) smaller width set by the value 
of J. This result is in good accord with the dispersion and width found for the low energy 
"foot" in recent angle-resolved photoemission (ARPES) data and t-J model results (there 
t 0.4 eV) []TU|. This is illustrated in Fig. 1(b), where the low-energy peaks in A{k,uj) are 
fitted by (full and dotted lines) E^ = A + J /2{coskx + cosky)"^, with A = 2.4t, rather than 
by the SDW (strong-coupling) result, i.e. Ej: = y^e^ + A^ = A + J^cosk^ + cosky)^. The 
overall agreement between the ARPES width and the Hubbard model data (for t ~ leV) is 
significant because it shows in fact that the energy scale of the low-lying insulating band is 
controlled by many-body effects beyond the mean-field SDW result. 

Our findings for the spectral-weight dispersion uj{k) are schematically summarized in 
Fig. 2(a), which emphasizes again the different energy scales (J, U, SDW ||2^) involved. 
These energy scales become particularly pronounced in the {U = 12t)-case in Fig. 3. 

Before moving to the doped situation, we note that the jSt = 3, U = 12t result in Fig. 3 
does have the valence-band maximum at the M-point and not at the X- or (7r/2, 7r/2)-points 
in contrast to the l3t = 10 results for both U = 8t and U = 12t (not shown). This at first 
puzzling "high-temperature" result reflects the fact that ai (3t = 3 the spin-spin correlation 
length ^ is about a factor of 2.5 smaller than the QMC lattice extension. The system then 
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shows precursor effects of a metal, which move the spectral weight (valence-band) maximum 
- in agreement with the metallic situation in Fig. 4 - to the M-point. Otherwise, the low- 
energy band is found to be essentially unaffected by changing temperature from T = O.lt to 
T = 0.33t. 

Keeping this in mind, we consider in Fig. 4 the low-energy electronic structure in the 
metallic regime for doping < n >=0.95. 4(a) shows again the 3-D plot of A{k,uj) and 4(b) 
the dispersion relation with the degree of shading representing the intensity of Aik, cu), as in 
Fig. 1(b). Like in the half-filled case, we observe two general features, which are both seen 
in recent photoemission experiments 0]: a broad "background" of ~ 4t — 6t spanning the 
lower and upper Hubbard band and a pronounced low-energy "foot" of significantly smaller 
width, which is clearly resolved between F^X and F— i'(7r/2, 7r/2). The situation is depicted 
schematically in Fig. 2b. Width and dispersion of the low-energy "band" are also in good 
accord with exact diagonalization results of 4 x 4 clusters |[1^,|15|. 

The results for the doped case have several important implications: First, they reveal that 
the lowest energy "band" in the insulator and the "band" that crosses the chemical potential, 
C(j = /i, in the hole-doped metal are rather similar: the low-energy band is separated from the 
broad valence-band background ( "LHB" (lower Hubbard band) in the schematic drawings of 
Fig. 2) in the same way; it has similar dispersion and it has a similar intensity modulation as 
a function of energy. Thus, the lowest-energy metallic band appears to be effected by similar 
many-body physics as in the insulating regime, namely magnetic correlations connected with 
the (now) short-range AF order. This is not in contradiction but instead substantiated by 
the fact that the metallic band develops its maximum at the M-point: as pointed out above, 
this happens as soon as the spin-spin correlation length ^ is smaller than the length (L=8), 
a situation obtained for the filling < n >=0.95. 

Furthermore, our QMC data display similarities to results derived by Kampf and Schri- 
effer 0. There a phenomenological model for the spin susceptibility was used to extract 
the leading-order contribution to the self-energy from AF-fiuctuations. The spin-spin cor- 
relation length ^ was an input parameter to the model and by varying this correlation 
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length the model evolved continuously from the insulating AF to the Fermi-liquid regime. 
For decreasing correlation length weight was transferred into the QP peak from incoherent 
backgrounds, which (for C, large) developed into the upper and lower Hubbard bands (our 
structures "UHB" and "LHB"). This is schematically summarized in Fig. 2(b). 

Another noteworthy feature of the doped, metallic situation is that the intensity change 
(but not the width) as a function of k for the higher-energy background in A{k,uj) still 
follows essentially the AF SDW picture with, in particular "shadow bands" resulting from 
the AF short-range order being clearly visible at T-^u — ~ — 5t) and M-{uj — ~ +7t) 
points. Even remnants can be detected of folded back "shadow bands" near the M-point for 
u < fi, which in the SDW-picture have much less oscillator strength and spectral intensity 
than the original band (between F— >X) [^. The findings confirm again to a certain extent 
the phenomenological work by Kampf and Schrieffer [^]. Finally, we would like to mention 
that an important detail of the QMC data, the rather extended flatness of the energy band 
near the X(7r, 0)-point, which is in good agreement with ARPES experiments of Dessau et 
al. for Bi 2212, has previously been resolved in QMC work both for the one-band |Q and 
the three-band [ PT| , P7| Hubbard models. This rather extended flatness, extending like the 
ARPES data not only into X— >r, but also into the X^M regions (as displayed in Fig. 3) is 
already inherent in the undoped low-energy structure near X, a fact which has recently also 



been found in 2-D t-J model studies by Dagotto et al. |jTO|. Its extension, in particular into 
X— >M direction, is not consistent with available one-electron band calculations It can 
be explained by a conventional self-energy diagrammatic analysis summing over the leading 



spin-fluctuation diagrams . It is thus a many-body effect related to magnetic correlations 



consistent with the arguments given in this work. 

In summary, we have studied the evolution of the 2-D Hubbard model from insulator 
to metal in terms of the electronic spectral weight, obtained from the maximum-entropy 
analytic continuation of QMC data. These results, combined with recent ARPES data, can 
be taken as strong indication that the QP dispersion of the high-Tc compounds, not only in 
the insulating limit but - particularly - in the metallic situation, has a many-body origin: 



the coupling of the quasiparticles to antiferromagnetic correlations. 

We would hke to thank R. Laughlin, S. Maekawa, A. Muramatsu, D. Poilblanc, H. Schulz, 
Z.-X. Shen and, particularly, D.J. Scalapino for instructive discussions. The calculations 
were performed at the HLRZ Jiilich and at the LRZ Miinchen. 
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FIGURES 

FIG. 1. Single-paxticle excitation for the 8x8 Hubbard model at half-filling, U = 8t, f5t = 10: 

(a) The single partiele spectral weight A{k,uj) versus lo and k; (b) u versus k "bandstructure"' , 
where sizable structure in A{k, u) is represented by strongly shaded areas and peaks by error bars. 

FIG. 2. Schematic plot of the bandstructure: (a) the insulator and the different energy scales 
(see text) involved; (b) for the metallic situation. 

FIG. 3. Single "bandstructure" u> versus k for U = 12t, f3t = 3. Again strongly shaded areas 
correspond to maxima of A{k,Lj). 

FIG. 4. The same as Fig. 1, but now for doping < n >=0.95 and /3t = 3, U = 8t: (a) A{k,Lo); 

(b) uj versus k "bandstructure" . 
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